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ABSTRACT 

Context. After the occurrence of the type cllb SN 201 ldh in the nearby spiral galaxy M 51 numerous observations were performed 
with different telescopes in various bands ranging from radio to y-rays. 

Aims. We analysed the XMM-Newton and Swift observations taken 3 to 30 days after the SN explosion to study the X-ray spectrum 
of SN 201 ldh. 

Methods. We extracted spectra from the XMM-Newton observations, which took place ~7 and 1 1 days after the SN. In addition, we 
created integrated Swift/XRT spectra of 3 to 10 days and 1 1 to 30 days. 

Results. The spectra are well fitted with a power-law spectrum absorbed with Galactic foreground absorption. In addition, we find 
a harder spectral component in the first XMM-Newton spectrum taken at t ~ 7 d. This component is also detected in the first Swift 
spectrum of t = 3 — 10 d. 

Conclusions. While the persistent power-law component can be explained as inverse Compton emission from radio synchrotron 
emitting electrons, the harder component is most likely bremsstrahlung emission from the shocked stellar wind. Therefore, the harder 
X-ray emission that fades away after t ~ 10 d can be interpreted as emission from the shocked circumstellar wind of SN 201 ldh. 
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1. Introduction 

Massive stars at the end of their lives have undergone phases of 
more or less strong mass loss. The evolution of the stars and thus 
also the way their supernovae (SNe) evolve, depend strongly 
on their mass loss rates. The interaction of the SN shock wave 
with the stellar wind creates a hot region around the SN site, in 
which electrons are accelerated and produce radio synchrotron 
emission. The hot plasma around the SNe can be observed in 
X-rays. Strong X-ray emission is in particular expected when 
the SN shock wave breaks out o f the star and starts propagating 
into the circumstellar material dChevalier & Irwir3l201 ll I2012L 
and references therein). Owing to prompt observations with the 
Swift or the Chandra telescope, X-rays from the interaction of 
the SN shock with the circumstellar matter have been detected 
for a number of SNe (e.g. SN 2006jc or SN 2010il.llmmler et al l 
2008l: IChandra et al.l2012l respectively) . ISoderberg et al.l d2008^ 
detected a transient X-ray source with Swift, which was then 
identified as SN 2008D in the galaxy NGC 2770. The X-ray out- 
burst was ascribed to the shock break-out of the SN. 

On May 31, 2011 , a supernova explosion was observed 
dSilverman et al.ll201 ll) i n the nearby galax y M51 located at a 
distance of 8.4 ±0.7 Mpc dVinko et alJl2012l) . It was classified as 
a type lib SN (lArcavi et al.ll201 lal) . Several radio, optical, and 
X-ray observations followed to study the SN evolution, which 
also allowed to constrain the nature of the progenitor suggesting 
a compact progenitor star with a radius of ~ 10" cm and thus 



* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 



the classification of SN 201 l dh as a type cllb SN dArcavi et al.l 
1201 lbt ISoderberg et alj|2012l) . 

ISoderberg et all d2012l) studied SN 201 ldh using Swift 
and Chandra data in X-rays and the Submillimeter Array, the 
Combined Array for Research in Millimeter-wave Astronomy, 
and the Expanded Very Large Array in radio in the following 
weeks. They showed that the observed radio emission is syn- 
chrotron radiation of electrons, accelerated in the forward shock 
of the SN explosion, while inverse Compton (IC) scattering of 
these electrons produced X-rays. They estimated that the break- 
out of the shock out of the compact progenitor must have oc- 
cured at « 4 x 10 11 cm with a rise time for the break-out 
pulse of fbr ~ 1 min. However, the shock break-out pulse was 
not detected by any X-ray or y-ray observatory. 

We report the detection of a hard X-ray component in the 
early spectra of SN 201 ldh observed with XMM-Newton and 
Swift. The comparison of the spectra taken at different times af- 
ter the SN event ranging from t — 3 to 30 d allow us to consider 
this component as emission from the interaction of the SN shock 
with the circumstellar material in the immediate surroundings of 
the progenitor star. 

2. Data 

2.1. XMM-Newton data 

Shortly after SN 201 ldh occured, two observations with XMM- 
Newton were initiated. The first observation with the ObsID 
0677980701 was performed at t « 7 d after the SN explosion 
from 201 1/06/07, 5:19:56 to 201 1/06/07, 8:30:40 (UTC) with an 
exposur e of ~ 1 1 ks. All European Photon I maging Cameras 
(EPICs, IStriideretall200U iTurner et al.ll200ll) were operated in 
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Fig.l. XMM-Newton EPIC-pn images ~1 d (left) and -11 d 
{right) after the SN explosion in the colours red (0.3 - 1.0 keV), 
green (1 .0 - 2.0 keV), and blue (2.0 - 8.0 keV). 



full-frame mode and used with the thin filter. The next observa- 
tion took place from 201 1/06/1 1, 5:05:03 to 201 1/06/1 1, 8: 15:38 
(UTC) at t sj 11 d with ~ 10 ks exposure. The EPICs were used 
in the same configuration as in the observations before. After 
filtering good time intervals, only 4 and 2.6 ks of useful data 
remained for the observations 0677980701 and 0677980801, re- 
spectively. 

Due to low statistics, the MOS 1/2 data are not useful for fur- 
ther analysis. Therefore, we created images from the EPIC-pn 
data in the bands 0.3 - 1 .0 keV, 1 .0 - 2.0 keV, and 2.0 - 8.0 keV 
(Fig.[TJ. As can be seen in the EPIC-pn images, the X-ray source 
at the position of SN 201 ldh is brighter on 201 1/06/07 than on 
2011/06/11 and also appears more blueish. The red-orange ap- 
pearing source in the center is the nucleus of M 51 together 
with an unresolved ultra-luminous X-ray source (ULX). The 
other brighter blue or green sources are also ULXs in M 51 
dDewangan et alj|2005l) . 

We extracted EPIC-pn spectra of SN 201 ldh in a circular re- 
gion with a radius of 20" around the optical position. The back- 
ground spectrum was extracted in a 30" radius circle close to the 
source, where no X-ray source was detected. Hereafter, we call 
the EPIC-pn spectrum of 201 1/06/07 spectrum XMM1 and that 
of 201 1/06/1 1 spectrum XMM2. 

2.2. Swift data 

SN 201 ldh was also observed with the Swift satellite in a 
large number of observations t o stud y the evolution of the X- 
ray emission. [Soderber g et alj d2012l) derived the X-ray fluxes 
of SN 201 ldh from Swift observations taken with the X-ray 
telescope (XRT) from 2011/06/03 to 29 and analysed the in- 
tegrated spectrum from the observations of 2011/06/03 to 17. 
Since we would like to study the Swift spectra at times compa- 
rable to the XMM-Newton observations, we collected XRT ob- 
servations of SN 201 ldh starting ~ 3 d after the SN explosion, in 
the time period from 2011/06/03 12:04:02 (UT) to 2011/06/10 
17:08:32 (total exposure time of ~ 44 ks) and in a later pe- 
riod from 2011/06/11 04:55:01 to 2011/06/30 12:33:50 (total 
exposure time of ~ 85 ks). The XRT data obtained in photon- 
counting (PC) mode were proc essed with the stand ard proce- 
dures (XRTPIPELINE v.0.12.6. lBurrows et al J 120051) . Standard 
grade filtering (0 - 12) and screening criteria were applied. 

Events for the spectral analysis were accumulated within the 
same circular regions as for the XMM-Newton spectra, i.e., with 
20" radius centered on the optical position of the SN for the 
source spectrum and from a source-free circular region of radius 
30" close to the SN for the background. We used version v.013 
of the response matrices in HEASARC calibration database 



(CALDB) and the corresponding ancillary response files created 
using the task XRTMKARF. We call the integrated spectrum of 
201 1/06/03 to 10 spectrum Swiftl and that of 201 1/06/1 1 to 30 
spectrum Swift2. 

3. Spectral analysis 

The spectrum Swiftl is merged from the data of 2011/06/03 
to 10 and thus represents an average spectrum of before 
2011/06/10. Similarly, spectrum Swift2 is an average spectrum 
of after 2011/06/11. Therefore, the four analysed spectra have 
the following chronological order: Swiftl, XMM1, XMM2, and 
Swift2, with Swiftl and XMM1 corresponding to spectra taken 
at similar times. 

3.1. The longer persistent soft X-ray emission 

We analysed the two XMM-Newton and two Swift spectra si- 
multaneously to search for changes in the spectral components. 
Pooled (l201ll) reported that the Chandra spectrum taken on 
201 1/06/12 can be fitted with a power-law spectrum with a pho- 
ton index of T = 1 .4 + 0.3, absorbed by the Galactic foreground 
Nk = 1.8 x 10 20 cm' 2 . For the Swift spectra from 201 1/06/03 to 
17,[Spderberg et al.l d2012l) determined a photon index of T = 0.9 
-1.8. We therefore fitted the spectra first with a single power-law 
model, also assuming foreground absorption by Galactic A^h = 
1.8 x 10 20 cm' 2 . The first two spectra Swiftl and XMM1 can be 
fitted with a lower photon index of F — 1.1 (1.0 - 1. 30 with red. 
X 1 — 1.1 and d.o.f. = 33, whereas the spectra XMM2 and Swift2 
taken after 2011/06/11 are fitted well with F = 1.8 (1.5 - 2.0) 
with red. x 1 = 1.1 and d.o.f. = 23. Therefore, the spectra taken 
on the first ~10 days after the SN is significantly different than 
the spectra thereafter. 

We also fitted all spectra with a single thermal 
bremsstrahlung model. While the temperatures fitted for the first 
two spectra (f < 10 d) are unconstrained (£7\, remSi swifti > 9 keV 
and fc7b rems ,xMMi > 22 keV, respectively), the two later 
spectra (t > 11 d) are both fitted well with temperatures of 
Wbrems,XMM2 = 3(1-18) keV and £r bre ms,Swift2 = 3(2 - 8) keV 
(red.^- 2 = 1.1 at d.o.f. = 56). 

3.2. The fading hard X-ray emission 

A likely origin of the additional harder component in the ear- 
lier spect ra is free-free emission from the shocked circumstel- 
lar wind dChevalier & Franssonl 12003). If we keep the power- 
law component and include an additional free-free emission 
component for the earlier spectra Swiftl and XMM1, the pho- 
ton indices of their power-law component become higher with 
F = 2.0 (1.6 - 2.5) (red. X 2 = U at d.o.f. = 32). Moreover, 
we obtain an additional absorbing column density of N\i ma 
= 7.2 (1.3 - 15.0) xlO 20 cm' 2 . The flux of the power-law 
component is the same for all four spectra with the unab- 
sorbed flux being F pow (0.3 - 8.0 keV) = 1.1(0.8 - 1.5) x 
10~ 13 erg cm -2 s' 1 . The temperature of the bremsstrahlung com- 
ponent is not well constrained and we can only determine a 
lower limit of kTb iems > 42 keV. Figure [2] shows the confi- 
dence contours for the parameters fe7\, rems and F. The unab- 
sorbed flux of the bremsstrahlung component is i*brems(0.3 _ 
8.0 keV) = 1.3(0.9 - 1.7) x 10' 13 erg cm' 2 s' 1 and 1.1(0.7 - 
1.6)xl0' 13 ergcm' 2 s' 1 forSwiftl and XMM1, i.e., for t < lOd 



1 All errors in this paper given in brackets are 90% confidence ranges. 
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Fig. 2. Confidence level contours (63, 90, and 99%) for the pa- 
rameters fclbrems an d r. 



and t — 1 d, respectively. While this component has disap- 
peared at t — lid with an upper limit of the unabsorbed flux 
of ^brems(0.3 - 8.0 keV) < 0.6 x 10~ 13 erg cirr 2 s _1 , the power- 
law component seems to represent emission that is persistent up 
to 30 d or longer. Figure [3] shows the first two spectra Swift 1 
and XMM1 with a model consisting of a bremsstrahlung and 
a power-law component (upper diagram), while the later spec- 
tra XMM2 and Swift2 are shown with the single power-law fits 
(lo wer diagram). 

Hmmler & L ewinl d2003h on the other hand discuss the pres- 
ence of a blackbody continuum of the shocked hot circumstel- 
lar gas. If we fit the spectra Swiftl and XMM1 with a black- 
body model for the additional hard component we obtain a tem- 
perature of IcTbb = 1.6(1.7 - 2.8) keV. The luminosities are 
Lbb(< 10 d) = 1.6(0.8 - 3.7) x 10 39 erg s" 1 and L BB (7 d) = 



1.6(0.8-3. 1)xl0 jy ergs" 1 for Swiftl and XMM1, respectively. 
From the fits of the spectra XMM2 and Swift2, we obtain the 
following upper limits: L BB (11 d) < 0.8 x 10 39 erg s _1 and for 
XMM1, L BB (11 - 30 d) < 0.3 x 10 39 erg s~\ 



4. Discussion 

Many SNe show a softening o f their X-ray spectrum days to 
months after the SN event fe.g Jlmmler & K untz 2005, and ref- 
erences therein). If the X-ray emiss ion is mainly due to IC pro- 
cess as proposed for SN 201 ldh by Soderberg et"ail d2012l) . one 
expects dLx/dE oc £'-(p- 1 )/ 2 ) with p being the spectral index 
of the energy of the injected ele ctrons responsible al so for the 
radio synchrotron emission (e.g. JChevalier et al]|2006l) . This, of 
course, is based on a simple assumption that the particle spec- 
trum does not change over a broad spectral range. In reality, 
some acceleration processes or particle losses can modify the 
spectrum. As the power-law fits of XMM1 and Swiftl vs. Swift2 
and XMM2 have shown the photon indices are significantly dif- 
ferent (r = 1.1 [1.0 - 1.3] and 1.8 [1.5 - 2.0], respectively) and 
indicate temporal change in the X-ray spectrum. Assuming that 
the particle spectrum does not change its shape significantly, the 
softening of the X-ray spectrum is most likely not caused by the 
change of the slope of the non-thermal spectrum, but rather by 
the change of the emission components. The X-ray emitting gas 
behind the blast wave of the SN shock is expected to have high 
temperatures of 100 keV or higher, while the gas behind the 
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Fig. 3. X-ray spectra of SN 201 ldh: stacked Swift spectrum for 
t < 10 d (Swiftl, black in the upper diagram), XMM-Newton 
EPIC-pn spectrum taken at t » 7 d (XMM1, red in the upper 
diagram), XMM-Newton EPIC-pn spectrum taken at t « lid 
(XMM2, black in the lower diagram), and stacked Swift spec- 
trum for 11 < t < 30 d (Swift2, red in the lower diagram). The 
spectra in the upper diagram are fitted with a combined power- 
law + bremsstrahlung spectrum, while the spectra in the lower 
diagram are fitted with a power-law spectrum only. The Swift 
spectra appear at lower countrate per energy bin due to lower 
effective area of the detector. 



reverse shock is cooler (-1-10 keV lImmlerll2003l) . Hence a 
likely explanation for the softening of the X-ray spectrum is that 
the reverse shock emission becomes more dominant. 

The analysis of the XMM-Newton EPIC-pn and Swift XRT 
spectra of SN 201 ldh have shown that the X-ray spectrum 
changes during ~10 d after the SN event. If we assume that the 
longer persistent X-ray emission component can be interpreted 
as thermal bremsstrahlung, we obtain a temperature of -3 keV 
after t *1 1 d, i.e., ~ 3 x 10 7 K, corresponding to a shock velocity 
of ~ 1500 km s _1 , which is far lower than what was measured in 
the op tical spectra at these times (~10000 km s -1 . lArcavi et al.l 
1201 lbl) . Therefore, this component is better to be identified as 
reverse shock emission or IC emission. In this case, the addi- 
tional hot component observed only until ~10 d after the SN 
event might be emission from circumstellar gas shocked by the 
forward shock. 
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A similar X-ray emission was also detected for SN 1993 J 
with ROSAT few days after the SN explos ion, which occured 
in the galaxy M 81 at a distance of 3.6 kpc (Zi mmermann et all 
Il994l) . This emission was fitted well with a power-law spectrum 
with F = 1.0+0.25 or a bremsstrahlung spectrum with kTb iems > 
7 keV, and decreased exponentially during the ~40 d in which it 
was observed. 

Assuming free-free emission for the hard component, we 
obtain a lower limit for the temperature of kT\, lems > 42 keV. 
The flux at t = 7 d is F bKms (Q3 - 8.0 keV) = (1.1 + 0.5) X 



10 13 erg cm 



i~ 2 s _1 , corresponding to Lb re ms(0.3 - 8.0 keV) » 
9.3 x 10 38 erg s _1 . Using the lower limit for the temperature, we 
obtain a lower limit for the shock velocity of 



l6kT bK 



3/77 



> 9800 km s" 



(1) 



with a mean mass per free particle for fully ionised plasma of 
m = 0.61 m v . Soderberg et al.ld2012l) estimated v s * 0.1 c, which 
is about three times higher than this lower limit. The XMM- 
Newton and Swift spectra that only extend up to ~10 keV do not 
allow to constrain the temper ature for a brems s trahlu ng model. 
Using v s » 0.1 c derived by Soderberg et afl (1201 2l) . we esti- 
mate that at t = 7 d, the shock must have reached a radius of 
R s « v s t * 1.8 X 10 15 cm. Therefore, from the bremsstrahlung 
spectrum at t — 7 d, we derive a down-stream density of n = 6.4 
xlO 7 cm" 3 at/? s « 1.8 x 10 15 cm. 

On the other hand, we can also estimate the stellar wind 
density based on assumption of the progenitor star. SN 201 ldh 
was classified as a type lib SN with a compact progenitor, 
making a Wolf-Rayet (WR) star a likely progenitor candidate 
dChevalier & Soderberg|l2010h . The number density of the wind 
at a distance r from the WR star is given by the continuity equa- 
tion: 



n(r) = 



M 



4njjm p r 2 v(r) 



(2) 



where M is the mass-loss rate, fi is the mean atomic weight of 
the w ind material, and v(r) is the /3-velocity law (ICastor et al.l 
1975): 



v(r) = Vc» 1 



bR* 



(3) 



R* is the stellar radius, Voo the terminal velocity, b = 0.9983 a 
para meter that ensur es that v(R t ) > 0, and /3 = 4 — 8 for WR 
stars dSchmutdll997l) . The volume integral of equation (|2]l over 
a shell with radii R\, r and R s centered at the position of the WR 
star gives the total number of particles: 



f n{ 

Jv 



N = I n(r)dr. 



(4) 



Thus, from equation |@}, it is possible to obtain the average num- 
ber density in that shell: 



n = N x 



3*(*s-0 



(5) 



ISoderberg et al.l d2012l) estimated a stellar radius of R* = 
10 11 cm, a shock break-out radius of R\, r = 4 x 10 11 cm, and a 
mass loss rate of M = 3 x 10~ 5 Moyr _1 . With a typical wind 
velocity of = 10 3 kms -1 and R s = 1.8 x 10 15 cm at t = 7 d, 
we obtain a mean stellar wind density of h « 8 x 10 5 cm -3 from 
equation ([5]). This calculation is based on the assumption, that 



the stellar wind density is highest close to the stellar surface and 
decreases with ~ r 1 . However, the stellar wind material might 
as well form a dense shell around the star at a certain distance , 
as was most likely the case for SN 2006jc dlmmleret al.l l2008). 
This SN showed an increase in X-rays about 100 days after the 
explosion which was interpreted as emission from shock heated 
shell at R s « 10 16 cm with a thickness of AR a 2 x 10 15 cm. A 
similar condition might have been the case also for SN 201 ldh, 
however, with a shell being located closer to the star, which re- 
sulted in an earlier rise and decay of the X-ray emission. 



5. Summary 

Two XMM-Newton observations of SN 201 ldh performed ~7 
and 1 1 days after the discovery revealed that the X-ray spec- 
trum changes significantly at about 10 days after the explosion. 
The analysis of the XMM-Newton EPIC-pn spectra and addi- 
tional stacked Swift/XRT spectra extracted from many obser- 
vations of similar periods supports the existence of two spec- 
tral components with the harder component disappearing after 
~10 days. The softer compo nent can be identifie d either as IC 
emission as suggested by Soderb erg et al.l d2012l) or as reverse 
shock emission with an unabsorbed flux of F so f t (0.3 - 8.0 keV) * 
1 x 10~ 13 erg cirT 2 s -1 up to t a 30 d. The flux of the hard 
component if we assume that it is bremsstrahlung emission, is 
Fbr ems (0.3 - 8.0 keV) = (1.1 + 0.5) x 10~ 13 erg cirr 2 s" 1 at 
t — 7 d, and thus comparable to the soft component, and de- 
creases quickly thereafter. This indicates that this emission has 



its origin in the circumstellar matter that extends to ~ 10 cm 
and was heated by the SN shock. 
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